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Finite geometry is used to underpin operators acting in finite, d, dimensional Hilbert space. Quasi 
distribution and Radon transform underpinned with finite dual afiine plane geometry (DAPG) are 
defined in close analogy with the continuous (d oo) Hilbert space case. An essential role in these 
definitions play projectors of states of mutual unbiased bases (MUB) and their Wigner function-like 
mapping onto the generalized phase space that lines and points of DAPG constitutes. 



, , I. INTRODUCTION 

o : 

, In this article we consider, in a d-dimensional Hilbert space, the underpinning of information theoretic related 
^ ' operators with finite geometry. 

o : 

' The general idea of such a relation was studied intensely, [l|, . Our approach is focused on a particular branch of 
_ finite geometrical system, [7 9, 26]: dual affine plane geometry (DAPG) whose main features are outlined in section 
III. The operators that we discuss are the operators involved in "Mutual Unbiased Bases" (MUB) 0, BE^, HI)- The 
study is confined to dimensionality, d=p, a prime, ^ 2. We itemize the relevant definitions/properties concerning 
I— i. these operators in the next section, section II, in our attempt to have this paper self contained. 

' The general idea of geometrically underpinned Hilbert space is as follows. A Geometry is defined by interrelation 
among points and lines. Points are considered the fundamental entities and lines are aggregates of points. We 
specify the rules for the so called (finite) affine plane geometry (APG) and its dual (DAPG) in section HI. A specific 
arrangement of points and lines abiding by such rules form a realization of the geometry. The association of Hilbert 
' space operators with geometrical points and lines define a geometrical underpinning of the Hilbert space operators. 
Q^l Consistency of an underpinning is attested to by the consistency of the operators implied physical results/relation 
'— ' . and the interrelations dictated by the geometry. The present study relates to the so called state reconstruction 'eJ, [31 : 
^ ' the Hilbert space MUB state projectors allow the diagonal elements (of all MUB states) of the density operator to 
K-^ . define in full the density operator. We illustrate this with a study of Radon transform [g^, HI] which, hitherto, was 
QQ confined to studies in continuum. Consistency may here be checked by state reconstruction via the inverse transform. 

(N 

\^ • The focal issue in the present analysis is finite dimensional Radon transform. This transform, for the continuous 

\ case (d — )■ cx)) [13, [11], is widely used both in down to earth application (30l - l32j mathematical studies [H, [sj], and 

^-H ' state reconstruction analysis [ill , [isl , [ssl [36| . In the present work, dealing with finite dimensional Radon transform, 

T-H ' it is, much like in the continuous case, analyzed via quasi distribution, viz Wigner-like function defined here in the 

, finite phase space-like points and lines of DAPG. 
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The Radon transformation [la, [13, [1^ involves angular variables and thus its formulation in finite (Hilbert) space 
dimensionality is somewhat intricate. We overcome this by adopting a "physical" approach whose rationale, within 
the continuum, d oo case is as follows. The Wigner function, WA{q,p) maps an operator, A in Hilbert space 
onto a c- number function in phase space [13, [121 • When the operator is the state density operator, p, the resultant 
Wigner function, Wo(q, p), is quasi distribution, i.e. it enjoys many attributes of a distribution (it may, however, 
become negative) [ll[, [ij, [ij, [4l| . A particularly attractive attribute that underscores its role as quasi distribution 
is its marginal with respect to the position q. Thus integrating along a vertical line (the p coordinate) for a fixed 
q=x': J ^^|^(5(x' — q)Wp{q,p), gives the probability for the system in a state p (i.e. whose quasi distribution is 
Wp(g,p)) to have its position value x'. In other words this integral equals the expectation value of the projector 
la:')(x'l that projects the eigen state, \x'), of the operator x. The Wigner phase space mapping of this projector is 
W\x'><:x'\{<l,p) — 5{x' — q). This continues to hold when we replace the position operator by arbitrary MUB state 
projector, \x' ,9){6,x'\, with \x' ,9) the eigenfunction of Xg = xcosO +psin9, eigen value x' [O, [H, [i3l ■ In this case, 
since the Wigner function W^^i g-^^g .^.^ {q,p) = 5(x' — qcos9 — psin0), the integral is along the line y = —qsinO +pcosO. 
Thus here too the marginal relates to MUB projectors [l[, [lEIill and we have for the probability, P{x',9 : p), of the 
state p being found in x'9){9, x'\ {C — cos9, S — sin9), 

P{x',9:p)^tr{p\x',9){9,x'\)^ f ^S{x' ^ Cq ~ Sp)W,{q,p). (1) 



In this form we recognize the marginal probability of the MUB projector's expectation value as the Radon transform 
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of the quasi distribution, Wp{q,p). An explicit account is given in section II. We use these observations as our guide 
for the definition of the finite dimensional Radon transform. Indeed the inversion of the transform is, in effect, a 
state reconstruction: it gives the state's quasi distribution, in terms of the probabilities P{x' , : p). 

In the finite dimensional (d=p, a prime ^ 2) Hilbert space we define a Wigner function-like mappings of operators 
onto lines and points of DAPG. We show that our Wigner function-like mapping of the density operator, Vp{j), onto 
DAPG lines, Lj, - has all the quasi distribution attributes possessed by Wp{q,p) in the d — oo case. The marginals 
of Vp{j) that give the probabilities of the system being in some eigen function of the finite dimensional MUB sets 
are now given as summation are along the DAPG line, Lj these being assured via the function j, (defined in the 
text, a designates a DAPG point and j a line.) which corresponds to the delta function in the d oo. Moreover, the 
marginals here, much like the continuous case, are informationally complete i.e. they allow the reconstruction of the 
corresponding (quasi) distribution which, in turn, [ill [l6j determine the state. Thus our finite dimensional Radon 
transform is the marginals of the quasi distributions pertaining to the projectors of MUB states in complete analogy 
with the d — > oo case. 

The inversibility attribute of the transform allows the extraction of the system's state from lines' summations that 
is informationally complete. For clarity we give a brief explicit review of this for the continuous case in section II. 
The finite geometry with which we underpin the theory, viz dual affine plane geometry (DAPG) is outlined in section 
III. MUB is defined in section IV. Section V contains the basic underpinning theory. Here and in the succeeding 
section we list some useful formulae. The formulation of the finite dimension Radon transform with its explicit 
inversion as well as the definition our quasi distribution are contained in the succeeding section, section VII. Section 
VIII introduces underpinning of finite dimensional Hilbert space operators with affine plane geometry, APG. In the 
last section, section IX, we summarize and discuss the results. 



II. MUTUALLY UNBIASED BASES AND RADON TRANSFORM - THE d ^ oo CASE. 

The rationale of our analysis is the relation between informationally complete measurements which we enumerate 
via mutual unbiased bases (MUB) and the sought after state of the system. In this section we outline the (known 
e-g- [Hi) approach for the continuous case thereby clarifying, we hope, the (known, [lll[l3|) result that the state, 
actually its Wigner representative function, is determined from the diagonal elements of the density operator with 
respect to (all) the states in all the MUB bases. To this end we briefiy review the notion of MUB their measurements 
and the extraction thereof the Wigner function of the density operator. 

Mutual unbiased bases, MUB, in concept were introduced by Schwinger jlO;] in his studies of vectorial bases for 
Hilbert spaces that exhibit "maximal degree of incompatibility". The eigenvectors of x and p, \x) and \p) respectively 
are example of such bases. The information theoretical oriented appellation "mutual unbiased bases" were introduced 
by Wootters Two complete, orthonormal vectorial bases, i3i, B2, are said to be MUB if and only if {Bi 7^ B2) 

y\u), \v) e Bi, B2 resp., \{u\v)\ — const. (2) 

i.e. the absolute value of the scalar product of vectors from different bases is independent of the vectorial label 
within either basis. This implies that if a state vector is measured to be in one of the states, \u), of Bi it is 
equally likely to be in any of the states \v) of any other MUB, B2- (The value of the may depend on 

the bases, Bi,B2, which indeed is the case for the limit d — > 00, the continuous dimensionality.) MUB are found 
to be of interest in several fields. The ideas are useful in a variety of cryptographic protocols [l^ and signal analysis [2[ . 

We now outline, for the continuous Hilbert space dimensionality, some salient MUB features [l^ [l^, IH, l45j . 
Consider the so called quadrature [lll[i^[l3] operator Jtg and its eigen states state \x',9), 

Xe\x',9) =x'\x',e). (3) 

We now show that the states, \x' , 0) form a complete orthonormal basis for the (rigged) Hilbert space, with the vectors 
labeled by x' and the basis by 9, and bases with different 9, < 6* < tt, are MUB. We begin by noting, that, as can 
be easily checked [H ill [H 



Xe = U\9)xU{9) = Cx + Sp, 



(4) 
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where 

Uie)^e-^'^'^; a^^{x + zp), (5) 

and C — Cos9, S = SinO. Thus the solution to Eq(I3]) may be written in terms of the eigenvectors of the position 
operator [H, Hi, Sl-Iil , 

\x,e) = u\e)\x). (6) 

This defines our phase choice [l5|. We now use the well known result for the harmonic oscillator propagator, |40|. 
to get the X representation solution, 

{x'\x-e) = ^ ^-^{[x-+x'']C-2xx') (7) 

This form ascertains [IB] that lime_>.o(a;'|a;; 0) — 5{x — x'), and \mvg^T,/2{x'\x;9) — e^^ /V^- This phase differs 
from the more common one 0] and leads to an expression that is symmetric a: cc', a property that facilitates several 
calculations. We now verify, by direct calculation, that the bases {|a;; 9)} and {|x'; 9')} with 9^9' are MUB: 

\{x'-9'\x-9)\ = \{x'\U\9-9')\x)\ - , ^ (8) 
i\ , I , /I i\ I V !\ i\ ^2^\S{9,9')\ ^ ' 

S{6,9') = Sin{9 — 9'). Thus the number \{x';9'\x;9)\ is mdependent of the vectorial labels x, x'. We used the relation 
U{9')W{9) = W{9 - 9') and Eq®. 

The Wigner function that represents an arbitrary operator A in phase space [HI, [13, 0, [13, li^] is given by, 

WA{q,p) ^ [ dye'Py{q~y/2\A\q + y/2). (9) 



Thus the Wigner function maps an Hilbert space operator to a c-number function in phase space. An important 
attribute of this mapping is 

trAB^ J ^WA{q,p)WB{q,p). (10) 

It can be shown [ll|, [l^ that the Wigner function of the density operator is real though not non negative (in general) 
this, with EallOp. instigates its classification as a quasi distribution in phase space. This completes our brief review of 
MUB and quasi distributions for continuous Hilbert space |l9i]. Now we turn to the state reconstruction and Radon 
transform role for this case. 

The probability of obtaining x' upon measuring the quadrature operator Xg, Eq.(|31) for the state p is, [ill [isl [T6| 

p{9,x')^tr{p\x';9){9,x'\)^ J ^Wp{q,p)W^,, ,e^{q,p), (11) 

where we used Eg. pTI)) to write the RHS of the equation. Explicit evaluating Wxf,{q,p), via Eq.®, gives for p{9,x'), 

p{0,x')^ J ^Wp{q,p)S{x'-Cq-Sp). (12) 

This identifies 32] p{9, x') as the Radon transform of Wp{q,p). Thus the Radon transform [2^ [2^ TZf{L) of a function 
f(x,y) is the integral of the function over straight line, L. The line in Eq. p^ is a line in phase space whose given by 
x' = Cq + Sp. (The full transform requires the integration over all parallel lines i.e. over all values of x'). Solving for 
Wp{q,p) in terms of the observables p{9,x') [28j : 

Wp{q,p) ^n-^\nwp){q,p) = n-^ p{9,x'). (13) 

Here TZ represents Radon transform. Direct calculation gives [ll|> [13, [H, 



N 1 r f°° , , dp(x,9/dx 
Wp{q,p) = / d9r dx' ' ' 

TTJo J-oo 



x — qC — pS 



(14) 



Thus inverting the Radon transformation is state reconstruction as the state p is determined once its Wigner 
representative is known [Tlj . 

Note that p{9,x') plays a double role: (a) As a marginal it is required via Eq. pUl) . that the RHS of Eq. pTj) contain 
the Wigner function of the MUB state projector, \x' ,9){9,x'\. Alternatively, (b), as the Radon transform of Wp{q,p), 
this is consistent because the Wigner function of the projector is a delta function assuring, in the integration, that x' 
equals Cq+Sp in phase space. 
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III. FINITE GEOMETRY AND HILBERT SPACE OPERATORS 

We now briefly review the essential features of finite geometry required for our study [l], [13, [13, [11] . 
A finite plane geometry is a system possessing a finite number of points and lines. There are two kinds of finite 
plane geometry: affine and projective. We shall confine ourselves to affine plane geometry (APG) which is defined 
as follows. An APG is a non empty set whose elements are called points. These are grouped in subsets called lines 
subject to: 

1. Given any two distinct points there is exactly one line containing both. 

2. Given a line L and a point S not in L {S 3 L), there exists exactly one line L' containing S such that Lf]L' = 0. 
This is the parallel postulate. 

3. There are 3 points that are not coUinear. 

It can be shown [13, [IH] that for d — p™ (a power of prime) APG can be constructed (our study here is for d=p) and 
the following properties are, necessarily, built in: 

a. The number of points is d^; Sa, a = l,2,...d^ and the number of lines is d(d+l); Lj, j = l,2....d{d + 1). 

b. A pair of lines may have at most one point in common: Lj p| = A; A = 0, 1 forj ^ k. 

c. Each line is made of d points and each point is common to d+1 lines: Lj = IJ^ S^, Sa = 0^=1 ^'j- 

d. If a line Lj is parallel to the distinct lines Lk and Li then Lk \\ Li. The d^ points are grouped in sets of d parallel 
lines. There are d+1 such groupings. 

e. Each line in a set of parallel lines intersect each line of any other set: Ljf]Lk — 1; Lj jj" Lk- 
The above items will be referred to by APG (x), with x=a,b,c,d or e. 

The existence of APG implies [U, [13, [ll|the existence of its dual geometry DAPG wherein the points and lines 
are interchanged. Since we shall study extensivebly this, DAPG, we list the corresponding properties for it. We shall 
refer to these by DAPG(y): 

a. The number of lines is d^, Lj, j — 1, 2....d^. The number of points is d(d+l), Sa, a = 1,2, ...d{d + 1). 

b. A pair of points on a line determine a line uniquely. Two (distinct) lines share one and only one point. 

c. Each point is common to d lines. Each line contain d+1 points. 

d. The d(d+l) points may be grouped in sets of d points no two of a set share a line. Such a set is designated 
by a' € {aU Ma}, a' = l,2,...d. {Ma contain all the points not connected to a - they are not connected among 
themselves.) i.e. such a set contain d disjoined (among themselves) points. There are d+1 such sets: 

d{d+l) 

U 

Ra 

Ra Ra' 

e. Each point of a set of disjoint points is connected to every other point not in its set. 

A particular arrangement of lines and points that satisfies APG(x), x=a,b,c,d,e is referred to as a realization of 
APG. Similar prescription holds for DAPG. 

We now consider a particular realization of DAPG of dimensionality d — p,^ 2 which is the basis of our present 
study. We arrange the aggregate the d(d+l) points, a, in a d • (d + l)matrix like rectangular array of d rows and 
d+1 columns. Each column is made of a set of d points Ra = Ua'eauM "^"'i DAPG(d). We label the columns by 
b=-l,0,l,2,....,d-l and the rows by m=0,l,2...d-l.( Note that the first column label of -1 is for convenience and does 
not designate negative value of a number.) Thus a = m{h) designate a point by its row, m, and its column, b; when 
b is allowed to vary - it designate the point's row position in every column. We label the left most column by b=-l 
and with increasing values of b, the basis label, as we move to the right. Thus the right most column is b=d-l. We 
now assert that the d+1 points, mj(b),b — 0,1,2, ...d — 1, and mj{—l), that form the line j which contain the two 
(specific) points m(-l) and m(0) is given by (we forfeit the subscript j - it is implicit), 

™(^) = ^{c~ 1) +m{0), mod[d] b ^ -1, 
m(-l) = c/2. (15) 

The rational for this particular form is clarified in Section V. Thus a line j is parameterized fully by j = 
(m(— 1), m(0)). We now prove that the set j = l,2,3...d^ lines covered by Eg. dTSt with the points as defined above 



d 

= [J Ra; 

a=l 

= U Sa'-, 

a'eaUAIc 

= 0, a a' . 
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form a DAPG. 

1. Since each of the parameters, m(-l) and m(0), can have d values the number of lines (P; the number of points in a 
line is evidently d+1. DAPG(a). 

2. The linearity of the equation precludes having two points with a common value of b on the same line. Now consider 
two points on a given line, to(6i), m(&2); bi ^ 62. We have from Ea.([T5l). (6 ^ —1, bi ^ 62) 

m{bi) — -^(c — 1) + m(0), mod[d] 
bo 

m{b2) = y(c-l) + m(0), mod[d]. (16) 
These two equation determine uniquely {for d~p, prime) m(-l) and m(0). DAPG(b). 

For fixed point, m(b), c <^ m(0) i.e the number of free parameters is d (the number of points on a fixed column). 
Thus each point is common to d lines. That the line contain d+1 is obvious. DAPG(c). 

3. As is argued in 2 above no line contain two points in the same column (i.e. with equal b). Thus the d points, 
a, in a column form a set Ra = Uo'eauAf ^o", with trivially RaC\Ra' = 0, a 7^ a', and [J^i^^^ Sa = UQ=i-^cf 
DAPG(d). 

4. Consider two arbitrary points not in the same set, Ra defined above: m(&i), m{b2) {bi ^ 62)- The argument of 
2 above states that, for d=p, there is a unique solution for the two parameters that specify the line containing these 
points. DAPG(e). 

We illustrate the above for d=3, where we explicitly specify the points contained in the line j = (m(— 1) = 
(l,-l),m(0) = (2,0)) 

m\b -10 12 

• • • (0,2) 

1 (1,-1) • (1,1) • 

2 • (2,0) • 

For example the point m(l) is gotten from 

m(l) = 1(2 - 1) + 2 = 1 mod[3] m(l) = (1, 1). 



IV. FINITE DIMENSIONAL MUTUAL UNBIASED BASES, MUB, BRIEF REVIEW 



In a finite, d-dimensional, Hilbert space two complete, orthonormal vectorial bases, Bi, B2, are said to be MUB if 
and only if {Bi 7^ B2) 



V|m), \v) e 61, B2 resp., \{u\v)\ = 1/Vd. (17) 

The physical meaning of this is that knowledge that a system is in a particular state in one basis implies complete 
ignorance of its state in the other basis. 

Ivanovic (20j proved that there are at most d+1 MUB, pairwise, in a d-dimensional Hilbert space and gave an explicit 
formulae for the d+1 bases in the case of d=p (prime number). Wootters and Fields [8] constructed such d+1 bases 
for d = 2™ with m an integer. Variety of methods for construction of the d+1 bases for d = are now available 
0, Hll, |23j. Our present study is confined to d = p 7^ 2. 

We now give explicitly the MUB states in conjunction with the algebraically complete operators [l^ [l^ set: Z,X. 
Thus we label the d distinct states spanning the Hilbert space, termed the computational basis, by \n), n — 
0, l,..d-l;|n + d) = \n) 

Z\n) = X\n) = |n + 1), w = e^27r/d_ ^^g^ 

The d states in each of the d+1 MUB bases [l^|2l|are the states of computational basis and 

1 '''^ 

|m;6) = — ^a;^"("-i)-"'"|n); 6, m = 0, 1, ..d - 1. (19) 
Vd p 
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Here the d sets labeled by b are the bases and the m labels the states within a basis. We have [2l| 

XZ''\m;b) =uj"'\m-b). (20) 

For later reference we shall refer to the computational basis (CB) by b=:-l. Thus the above gives d+1 bases, b=- 
l,0,l,...d-l with the total number of states d(d+l) grouped in d+1 sets each of d states. We have of course, 

{m;b\m';b) ^dm,m'; |(m;6|m';6')l = b ^ b' . (21) 

This completes our discussion of MUB. 



V. GEOMETRIC UNDERPINNING OF MUB QUANTUM OPERATORS 

We now consider a DAPG as underpinning a two sets of operators acting in a d-dimensional Hilbert space, these 
are 

i„; a = l,2....d(d+ 1), Pf, j^l,2,..(f. 

Here Aa are associated with the d(d+l) points, Sa while Pj are associated with the (f lines, Lj. 
We now define interrelation among the operators of the sets in a way that reflects the geometry, 

^ = ^E^.- (22) 



d 



These entail, using DAPG(a,b,d), 



T.^<^~lY.Pr- (23) 



d 

aej j 



We now list some known finite dimensional features. These will allow its underpinning with DAPG which are 
presented in the succeeding Section. 

We consider d=p, a prime. For d=p we may construct d+1 MUB [H, 0, HO) HH • Returning to our state labeling we 
label the MUB states by |m, b). We designate the computational basis, CB, by b=-l, while 5 = 0, 1, 2...d— 1 labels the 
eigenfunction of, resp. XZ''. m labels the state within a basis. (Note that assigning b=-l to the CB is for reference 
only.) The projection operator defined by, 

Aa = \m,b){b,m\; a = {b,m}; & = -1, 0, 1, 2...d - 1; m = 1, 2, ..d. (24) 

The point label, a — (to, b) is now associated with the projection operator, Aa ■ We now consider a realization, 
possible for d=p, a prime , of a d dimensional DAPG, as points marked on a rectangular whose width ( x-axis) is 
made of d+1 column of points, each column is labelled by b, and its height (y axis) is made of d points each marked 
with m. The total number of points is d(d+l) - there are d points in each of the d+1 columns. We associate the 
d points m — 1,2, ...d, in each set labelled by b ^q,; a ^ {m^b) to the disjointed points of DAPG(d), viz. for fixed 
b a' e a U Ma form a column. The columns are arranged according to their basis label, b. The first being b=-l, 
a_i = (to, —1); TO = 0, 1, ...d — 1, signifies the computational basis (CB). The lines are now made of d+1 points each 
of different b. To a line Lj we associate an operator Pj. Now DAPG(c) (and Eq. ip^ ) implies: 

d d 

^\m,b){b,in\ = ^ Aa' =i 

m a'^a\jMa 
d(d+l) 

^ Aa ^ {d+l)i. (25) 
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Returning to Eq. ([23)) we have, 



Y^p, = di. 

3 

d+l 

= (26) 



These imply 



<r{i„P,} = A„,,; A„,, = r' " ^ (27) 

0, a 9 J. 



To prove Eq.(l27l) we note Eq. ([HI [Ml to write, 

■^aPj = Aa + ^ ' Aai Aq, — A^ 01 G j 





1 







1 


2 


ic/2-. 


=0, 


-1) 


^(0,0) 


Ao,i) 


^(0,2) 


ic/2-. 


= 1, 


-1) 






Al,2) 


(c/2-. 


=2,- 


-1) 


^(2,0) 


-4(2,1) 


^(2,2) 



= ^ ia'ia - « 9 j, (28) 

a' 

taking the trace imphes (j27p . This result, Eq. (j27p . involves geometrical factors only. We note that the Aaj may 
equally be viewed as 

fr{i„P,} = A„,,; Aa,, = j^' ^. ^ • (29) 

[0, J 9 a. 

The proof for this case is given in Appendix A., 

e.g. for d=3 the underpinning's schematics is (the choice of c/2 for the CB vectors will prove convenient below). 

m\b 

A 

1 A 



The geometrical line, Lj,j = (1,2) given above (end of Section III) upon being transcribed to its operator formula 
is via Eq. ipSl) , 

Pj{c/ 2=1, ma=2) = ^(c/2=l,-l) + ^(2,0) + ^(1,1) + ^(0,2) " I- (30) 

Returning to Eqs. (|24ll9p . these equations imply that, the projection operators A^, in the CB representation are 
given by, 

. . _/^; s = (n-n')(|[n + n'-l]-m), 5^-1, 

First we argue that every {Aa=m.b) in column b has in every other column b' 6 7^ 6', b,b' ^ —1 one pro- 
jector such that {Aa=m.b) = {Ac' =m' ,b') ^ ^, '■ Considcr two distinct columns, b,b' {b,b' ^ —1) and given the 
matrix elements n,n' {n 7^ n') of a projector {Aa=m,b) „/, compare it with [Aa' =m' ,b') If s fEq. pip i is 7^ s' i.e. 
|(n + n') — 1) — TO 7^ ^(n + n') — 1) — to' pick another projector in the same column, b' (i.e vary m'). Since m' = 
0,l,2...d-l there is one (and only one) {Aa' =m' ,b') ^ such that (Aa=m,b)„ — (^a'=m',&')n ^O'^ consider another 
matrix element (^Aa)^-,- We have trivially that [Aa)--, = (Aq,')_ iff n+n' = n + n' . i.e. all matrix elements (n,n') 
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with n+n'=c (constant) are such that (^q),^ ^, = (^a')n These elements are situated along a line perpendicular 
to the diagonal of the matrices. We refer to this perpendicular as FV (foliated vector), it is parameterized by c. 
We now assert that all other (non diagonal) matrix elements are unequal. i.e. for b ^ b' , (^Q=m,6)„„/ 

7^ {Aa'=m' ,b') n n" ^ 3 FV. Proof: Let two elements n,n' and 1,1' with n ^ n'; I ^ V in the two matrices 

be equal. Thus (c=n4-n', c'=l+r): 

b, ^ h' , , 

-(c— 1) — TO — — (c— 1)— TO, and 



= |(c'-1)-to'. 



(32) 



These imply c=c', QED. Now consider s=0. Then all the matrix elements along FV are 1/d. We have then that for 
(^")n n' ~ i^°')n n' ~ ^'^ /'^' ^ 7^ 0, d-1 matrix elements along FV are all distinct. The diagonal is common to all. 
We have, thus, a prescription for d projectors, ;,), one for each b, {b ^ —1), all having equal matrix elements 
along FV labelled by c. We supplement these with the projector A(c/2,-i) = |c/2)(c/2| to have the d+1 "points" 
constituting a line j . (|c/2) being a state in the CB.) Thus our line is formed as follows: It emerges from A(c/2,-i) 
continues to ^(m(o),o) in b=0 column. Then it continues to the points ^m(b),b in succession: b=l,2...d-l with 
m(b) determined by 

^(c - 1) - to(&) = ^(c - 1) - to(6 + 1). 

Thus the two parameters, c=2m(-l) and m(l), determine the line i.e. j=(m(-l),m(l). The general formula for the 
line is thus Ea. (jl5l now acquiring a meaning in terms of the point operators, ^Q=m(6).6- It is of interest that, if we 
associate the CB states with the position variable, q, of the continuous problem and its Fourier transform state, viz 
b=0 (cf. Ea. (jl9p ). with the momentum, p, we have that the line of the finite dimension problem is parameterized 
with "initial" values of "q" and "p" i.e. m(-l) and m(0). 

The discussion of the properties of the line thus defined confirm that these lines realize DAPG lines. The analysis 
above indicate that the line operator, Pj={m{-i),m{o))- We now list some important consequences of this. We have 
shown that the matrix elements along a FV direction are the same for all the point operators Aa^j. Indeed that is 
how we defined our lines. On the other hand we argued that the matrix elements not along the FV are all distinct. 
Thence summing up d such terms residing on a fixed Pj {excluding the b=-l and the diagonal term) sums up for each 
matrix element n,n' the d roots of unity for all matrix elements not on FV, hence for all c. 



: E 



Ac.-l)„,^,=0;n,n'3n + n' = c; a_i = |c/2)(c/2|. (33) 



Thus (-Pj)„ = ( 'J2aej ~ n' ^ "^^^ aloug FV, and is 1 along the diagonal at c/2=m(-l). The sum over the 

matrix elements on a FV, which are the same for all the Aa^j^^-i simply cancel the l/d. We illustrate this for the 
example considered above Ea.([5D|). viz: d=3, line with m(-l)=l, m(0)=2, i.e. j=(l,2) Evaluating the point operators, 

Aa , 






Ac/2=1-1) =\^^ly '^(2.0) ^ Q I ^„ 1 ^1 '^(1-1) = ^ I 1 M 'Ao,2) = ^ I 1„ \ (34) 

and evaluating the sum, Ea. (pO| . gives 

-Pj:(m(-l) = l,m(0)=2) = | ^„ ^ ^ | . (35) 

Quite generally, 

(P , . , _ /^"^"""'^"^°^'5{(„+„,),2™(-l)} . . 
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Thus, 

iP!=c/2,mo )",«' = Sn,n' . Z.e. ^ / Vj. (37) 

Theorem : trPjPji — dSj ji. 

Proof : Eq.^ ^ trPf = trl = d, 

Eq^ -> irPjPj = tr{Pj ^ i^) - trP,- 

aej' 

Eqm -> irP.i^e, - 1 = 0, / ^ j. 
Where we used trPj = 1 as follows from Eg. ([57)1 and that two distinct lines share one point (see DAPG(b)). 

VI. FURTHER ATTRIBUTES OF THE LINE OPERATOR, Pj 

The line operator, Pj is hcrmitian being the sum of hermitian operators. Noting that, trivially, trPj = 1 and 
P^ = I, Eq (P7|) implies the spectral representation 

^73^_^'^".»' n = l,2,...^ {^elements), 

\^3) — \ r „ d+3 d+5 J td-1 „7 y"^"! 



{^elements.) 



We thus have that 



is a projection operator onto the eigenstates of Pj with eigenvalue +1. Recalling that Pj — X^aej ^ ^ ^^"-^ 
P^ = I implies (what we term) the Fluctuation Distillation Formula (FDF): 

The proof is given in Appendix B. 

Summing both sides of the Eq. (P7|) over a and use Eg. ([25)1 and Eg. ipT]) to write, for the LHS 

= tr[{P,+i)P.j\=tr[P,+i]. (41) 

Where in the last step we used the relation, Eg . ([371 . Pj — I (valid for ^ = |c/2)(c/l|). 

To illustrate the spectral decomposition we consider again the line j — {to(— 1) = 1),to(0) = 2} in d=3. This line 
runs through the points (1,-1), (2,0), (1,1) and (0,2). The line operator -P(l,2) was given above, 

/O a;\ /1/2 uj/2\ / 1/2 lj/2\ 
Pi,2) = 010 = 10 - 10 . (42) 
0/ \ujy2 1/2/ \-ujy2 1/2/ 

Thus the matrix is diagonal in the orthonormal basis, 

0J/V2 \ /0\ / L0/V2 



^/V2) \0 J V-cjVV2 
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The associated projection operator, 

P(i,2) = = -^[c^lO) +c.|2)]i=[u;(0| +0.2(21] + |1)(1|. (43) 

Let Mj be a d(d+l) matrix. Its d+1 columns are made of the d elements, {n\m{b),b) with \in{b),b) the state 
whose projection, \m{b) , b) {b, m{b)\, is the point of the line j in the column b. Define Mj as the d columns matrix 
whose rows are made up of d+1 columns corresponding to the adjoint of Mj. For example, using Eq.(|30|) the line 
j : c/2 = 1, m(0) = 2 is given by the aggregate of points: (c/2 = 1, to(0) = 2), (2, 0), (1, 1), (0, 2) giving as a normalized 
state line operator. 



/O 1/V3 1/V3 1/V3 , . 
-j==\ 1 1/V3 ,M] = -=== 

V Loyvs LoyVi Loyvs ' 



/ 1 \ 

1/V3 wVv^ w/a/S 

1/V3 UJ/V3 uj/V3 

Vl/V3 1/V3 Lu/Vs) 



While the line j — c/2 = 0, mo = : viz 

(c/2 = 0, m(0) - 0), (0, 0), (1, 1), (2, 2) => M(,/2=o,m(o)=o) 
The construction of the matrices assures 



1 

7t 



1 1/V3 1/V3 1/V3 

1/^/3 wVVs w/Va 

1/^/3 wVV3 w/V3, 



trM]Mj 



1 



trM 



1 



d+1 



J + j'- 



(44) 



The proof is based on the DAPG attribute that two distinct lines share precisely one point, and for equal b distinct 
states are orthogonal, e.g. 

*^M[^/2=0,m(0)=0)^(c/2=l.m(0)=2) = \- 

The normalized "line" matrix, Mj, is the "square root" of the normalized "line" projection operator, viz: 

(45) 



' ^ d+l 
This is almost self evident: The RHS equals 

^A„=^|m(6),6)(6,m(6)| 

a^j b 

While the the LHS is (for the n,n' matrix element) 

Y,{n\m{b),b){b,m{b)\n'). 

b 

i.e. they are identical. 

We now demonstrate that these line operators are geometric in origin. Thus we associate the line operator Mj 
with Mj, defined by the replacement, in the former, of every column b (7^ —1), as follows. Instead of the elements 
{n\m{b),b), n = 0, 1, ...d — 1 with m(b) a point in the line j, we have in the column b of Mj, 1 at the row corresponding 
to m(b). Thus in our example of j=(l,2) for d=3 (cf. Eq. (pO|) ). we have. 



^ /O 1/V3 1/V3 l/VS \ ^ /O l\ 

M(i,2) = ^ 1 ^/V^ wVV3 l/VS ^M(i.2) = ^1010. 
^4 V LuyV3 ujyV3 0.2/73 J V4 \^ 1 y 

Mj is geometrical: it involves, in essence, the drawn line. We note that the transition between Mj and Mj is "local" 
unitary transformation: we require a distinct unitary transformation for each column (basis) ,b. 
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VII. DAPG UNDERPINNED QUASI-DISTRIBUTIONS 



We now define a Wigner function like quasi-distribution, Vp[j — (m(— 1), m(0))]. It maps the density operator p 
onto the lines, j = {(m(— 1), m(0))}, and points a = {m(6)}, of DAPG. This, Vp{j), Wigner function-like completely 
determine and is determined by the state of the system,/),. ( An essentially equivalent definition relates to arbitrary 
operators in the d- dimensional Hilbert space.) We then consider summation of Vp{j) along a fixed values of a that 
represent an MUB projector, as the finite dimensional Radon transform of Vp{j). The inversion of this, viz. the 
reconstruction of Vp{j) from summations along such points is our finite dimensional inversion of the Radon transform. 

Returning to the line operators, Pj, j = l,2,...cP, j = (m(— 1), to(0)), Eq. lfTS]) . we utilize their orthogonality, 
trPjPji = dSj_ji, to write, 

P= {^^PP^)Pj = E (46) 

j 3 

Where we defined Vp{j) = tr(^pPj). A partial list of attributes of Vp{j) is the following. 

i.Vp{j) = {Vp{j)y. 

2- \ Ef VpU) = trp{} Ef Pj) = 1, CfEq.m- 

3- trp = i Ef {trpP,)trP, = i Ef (trpP,) - 1- 

4. trAB = i Ef {trABP,)trP, = i Ef Vab{j). 

5. trAB = ^ Ef (trAP,) Ef {trBp,,)trP,P,, ^ \ Ylf ^A(j)T^B(i). 

In finite dimensional studies it is convenient to use unitary operators [lol . [l8| . Thus the probability, given that the 
system is in the state p, to measure it to be in the state \m, h) i.e. to be in an eigen function of XZ^ with eigenvalue 
without regard to any other probability is ir(p|m, 6) (5, to|) . Here |m, 6) is the eigen-function of the unitary operator 
XZ^, cf Eq. ipp]) . ir((o|m, 6) (6, rn|) is the finite dimensional observable that corresponds to tr{^p\x\6){9,x'\) of Eq. 
(|12p in the continuous case. Note that we may regard the c-number function gotten upon mapping, a la Wigner, the 
operator p\x' ^ 0) {9 , x' \ onto phase space as a marginal quasi distribution of PVp(<?,p) cf. Eq. (fTT|) . In what follows we 
introduce, for the finite dimensional case, Wigner like mapping, now onto DAPG coordinates, that, correspondingly, 
relates its marginals to the full quasi distribution for the MUB state projectors. 
Recalling, Ea. (|24l) . A^ = |m, 6)(m, 6|, a — {m,b), we may write, 

1 

trpA^^-Y,Vp{j)VAAj), cf. 4 above. (47) 



d 



Noting that. 



Va^ (j) = trAo^Pj = A„j cf. Eq.m- (48) 

i.e. 

trpA^ = ^''PPj = ^ E i^PP3^<^,3 - ^ E ^P^^)^A^ (j)- (49) 

jea j j 

These equations correspond to Eas. pT|) . (IT^ of Section II. Thus we identify trpA^ as the (finite dimensional) Radon 
transform of Vp{j) - it sums up the values of Vp{j) for j E a. (The Aaj plays the role of the delta function.) These 
equations correspond to Eq. (ITTl [TU)) . 

To obtain the (finite dimensional) inversion to the transform we consider, 

E Hp^'-) = ^ E ^"-pp^ E ^".^ ^ ^ E E ^"-pp^ = 1 + ^p(j)- (50) 

aej' j aej' j€a oiGj' 

Where in the last step we used the DAPG based relation, 

E E = E Pr + dPj + dp,. 

(aej){j'ea) j' = l 
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Thence, the inversion is 

^p(j) = E*KM")-i- (51) 

This could have been gotten directly from the definition of the operators however we could perhaps miss thereby 
some of the insight that the lengthy derivation provides which, in turn, underscores its relation with the continuous 
inverse Radon transform, Eq. (|13|) . 



VIII. AFFINE PLANE GEOMETRY (APG) FORMULATION 

The central work in the underpinning of finite dimensional MUB operators with finite geometry, T], is given in 
terms of lines and points of affine plane geometry (APG) rather than our choice of DAPG. An advantage of this, 
APG, scheme is its apparent similarity with the d —> cxd case in that it involves square arrays and states projectors are 
straight lines (albeit in a modular sense) in a "discrete" phase space. It allowed a direct imposition of translational 
invariance and extension to dimensionality d=p"; n > 1, p a prime, In this section we recast our DAPG 
underpinning into an APG one by interchanging lines and points and in particular the symmetrical meaning of A^j-, 
Eq. (I27|) . (|6ip is shown to allow the formulae for the Radon transform remain intact. 

Recall that within the DAPG underpinning a line was defined by the two points: m(-l) and m(0). The first, m(-l), 
was associated with modulated position, as it relates to the eigen values of Z, viz the computational basis states. We 
shall refer to it by ^. The second, m(0), was associated with (modulated) momenta - it being the eigen values of X 
(b=0, Eq. ipO)) . i.e. the Fourier transform of the former states. We shall refer to it by 77. Thus a DAPG lines are 
j = (m(— 1) = ^,m(0) = 77); ^,rj — 0,l,2,...(i— 1. Now consider a. d ■ d square array, d — p, aprime ^ 2, whose 
(discrete) coordinates along the x axis is labelled by ^ and the y axis by rj. We interpret each point, a, in the array, 
a = (^, 77) as underpinning a DAPG line j. Thus the "image" of each DAPG line is a APG point. We now consider 
lines in this array: 

r/ = r^ + s; mod[d\; r, s = 0, 1, 2, ...d — 1. 

C = s'; mod[d]; s' = 0, 1, 2, ...d - 1. (52) 

Eg. dS^ defines d(d+l) "straight" lines: there are d^ possibilities for r and s and d for s'. Each line contain d points: 
(^i:'/i): * = 0,l,...c?— 1. The line xi = s' contain d points as well: {s',r]i), i = 0,l,...d— 1. this proves APG(a), of 
section III. Similarly, the proofs of the validity of APG(x), x=b,c,d and e for the lines, Eq.(??) and points forming 
the array are trivial, e.g. to prove APG(b): Consider two distinct lines, rj — ri^ + si, rj — r2^ + S2, fi 7^ ''2, si 7^ •S2- 
Let this lines share a point {^o,Vo)- This implies, rifo + si = '"2'?o + S2- This implies a unique {£,Q,rjQ) : 770 = 
T^iCo + si; = (si — S2)/{r2 — ri). For ri = r2, si 7^ S2, no common point exist (the lines are parallel). Thus the 
square array with points labelled by (^, 77) and lines given by Eg ([5^ form a realization of APG. It is specified that a 
line j — {m{—l),m{0)) G DAPG is mapped to a point (^,77) € APG. (Note that a point, a, of DAPG underpins an 
MUB projector: a = \m,b){b,m\.) 

Theorem: In a DAPG realization, the d lines of DAPG that form the image of the d points of a APG line, 77 = + s, 
share a point. 

Proof: Consider an APG line, 77 = + s. It contain the d APG points {^i,rii — r^i + s), i ^ 0,1,. .d — 1. 
Now pick two arbitrary points i,i'. Their images in DAPG are the two lines (m{—l) = ^i,m{0) = r^i + s) and 
{m{—l) = ^i', 771,(0) = r^i' + s). The equation for the point they share, cf Eq. p^ . is 

^(2^, - 1) + r^, + s = ^{2^,, - 1) + r^,, + s mod [d], 

^ (b + rM^ - (^') = mod [d]. (53) 

this is independent of S,i i-G- the lines, i,i'=l,2,..d-l, share a common point at b=-r mod [d], thence the point 

is m(b=-r)=:r/2+s. For the APG line £, = s' the common point within DAPG is, trivially, at b=-l: m(-l)=s'. 
To illustrate the steps involved we consider d=3 with APG line given by 77 = ^ + 1. Thus the APG points involved 
are: (0,1), (1,2) and (2,0). Via Eq(|35]), Eq.® we have, 

'1 0' 

^^(0,1) = 1 0^ 'jj I ; Pii.2) = I 0. 1 1 ; ^(2,0) = | 1 1 . (54) 
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Now, Eq.dlSl) relates these to the DAPG point (i.e. the MUB projector) |0, 2)(2, 0| = A(o,2) 

1 





(55) 



where the last matrix is A(o,2)i cf. Eq. ([M)) 



To distinguish between DAPG and AFG underpinned operators we adopt the following scheme. The projector onto 
the MUB state \m,b)^ a point ala DAPG, was designated by ^(Q=m(&)) with the subscript indicating its coordinates. 
It is a line ala AFG, and will be designated by the subscript now gives the APG line's equation. The DAPG 

line operator Pj will, within APG underpinning scheme, be specified by the point P(^.ri) giving its coordinates. Within 
this notation our Wigner function- like mapping function, Ea. (j46p . is, 

Vp(f,77) -UVpQ(^,,). (56) 

Here the variables, ^,7y, signifies its APG underpinnings. In conformity with APG we have, (in correspondence with 
the DAPG relation, Eq.()) 



B 



For the sake of clarity we wish to iterate the relations between the Hilbert space operators Pj and Q^,ri'- Pj is a 
line operator within DAPG underpinning, j=(m(-l),m(0)) - the line is fully parameterized by the "position" variable 
value, m(-l) and the "momentum" m(0). Q^^n is a APG underpinned point operator with ^,7y the point's coordinates 
with ^ = to(— 1), r/ = m(0). The APG underpinning was constructed by identifying, 

Pj = {m{~l),m{0y) = Q(C=m(-l),r;=m(0))- 

Thus although the subscript of P designate a line underpinning (DAPG); while that of Q, is a point underpinning 
(APG), they are equal when the subscripts are numerically equal. 

The essential issue in the following is the simple, yet long winded, equivalence of our two accounts for the MUB 
operator, \m,b){b,m\: 

(a) As a DAPG underpinned point, a. In this case a — m{b), specifies the coordinate of the point in the d ■ (d + 1) 
points array. It is given by the point underpinned operator Aa. 

(b) As a APG underpinned line, A. In this case A = ?7(^) is the equation of the line whose constituent points, 
^,7] e A, correspond to the DAPG underpinned lines parameterized with m(-l), "position", and m(0), "momentum" 
coordinate: m(— 1) ^; m(0) rj. It is given by the line underpinned operator B^. 

The mapping of the line operator, Bx, is the Lambda function. 



(58) 



i.e. it is non-vanishing only with the point (^, rj) in the line A. It is a straight forward matter to show that the 
mapping of the density operator onto APG phase space like lines and points, Vp(^,77), is a quasi distribution - i.e. it 
possess the equivalent attributes, 1-5, heeded by our DAPG mapping given in the last section. AS an example we 
prove item 5: 



= ^ E (^^^Q^.-y) [trBQ^^,,) - i ^ Va{L vWeit v)- (59) 

The Radon transform of the quasi distribution, V^(^,?y), is trpB\ - since in the present case the line underpinned 
operator is the MUB state projector. Thus, 

^ d ^ 

trpBx = ^ E ^"^PQ^^v = ^E (^''^'55,»,)^(«,7?),A- (60) 



14 



This expression is in complete analogy with Radon transform in the continuum: we sum the quasi distribution over 
points 77) on the line A. Given that the system is the state p, the probability to measure it to be in |m, b), is given 
by Eg. ([60]) . This equation is analogous to Eq. ([TT|) : It expresses the probability of being in the MUB state (here \m, b)) 
in terms of summation of the quasi distribution (here trpQ^^rj) along a line determined by the mapping of the MUB 
projector onto the underpinning points (here A(^^ ,j) x = trB\Q^^ri)- 

IX. SUMMARY AND CONCLUDING REMARKS 

The finite dimensional, d, density operator, and mutual unbiased basis (MUB) states' projectors were mapped onto 
points and lines of finite geometry. These mappings were shown to be analogous to the Wigner function mapping of 
the density operator and MUB state projectors in the continuum, d ^ 00. In the latter {d — 00) case, the expectation 
values of the MUB state projectors expressed via the appropriate Wigner function scheme were observed to be the 
Radon transform of the Wigner function of the state itself. Inverse Radon transforms is, thus, state reconstruction in 
terms of the afore mentioned expectation values. 

The proposed finite dimensional map of the density operator possess the quasi distributional attributes over its 
underpinning geometry as the Wigner function over phase space. The mapping of the MUB state projectors are, like 
their Wigner function counters, lines. These lines are arranged as straight lines for the affine plane geometry (APG) 
underpinning. The underpinning with the dual affine plane geometry (DAPG) which was most extensively employed 
allows a simpler formulation. 

The expectation values of the MUB state projectors were used to define finite dimensional Radon transform. It involve 
summation along a line of the underpinning geometry. The inverse finite dimensional Radon transform is used for 
state reconstruction in close analogy with the continuum analysis. 

A brief summary of dual affine plane geometry (DAPG) in finite dimension, d, was given. The geometry was 
used to underpin projectors of states of mutual unbiased bases (MUB) scanning a d-dimensional Hilbert space. The 
dimensionality studied were d=p, a prime, ^ 2. 

The Wigner function, Wp{q,p), may be viewed as a mapping of the density operator that act in Hilbert space 
onto a c- number function in phase space. A finite dimensional, d, Wigner function- like, Vp{j), was defined a a 
mapping of the finite dimensional density operator, p, onto c-number function of lines, Lj = l,2,...d^ and points, 
a = l,2...d{d + 1), of dual affine plane geometry in d- dimensions. Vp{j) posses all the attributes of Wp{q,p) that 
qualify it as quasi distribution. A particularly attractive attribute of the Wigner function that underscores its role 
as quasi distribution is it marginals. In particular the marginals of mutual unbiased basis (MUB) state projector is 
itself a quasi distribution and is recognized as a Radon transform of Wp{q,p). This transform involves integration 
along a line in phase space which relates to the MUB state projector map onto this space. In close analogy the 
marginal of a finite dimensional MUB state projector involves summation along a DAPG line ans is itself a quasi 
distribution. This close analogy led, through its physical interpretation, to circumvention the angular involvement 
within Radon transformation and allowed a convenient definition of finite dimensional Radon transform concomitant 
with its inverse i.e. the state reconstruction. 

In closing we present in a comparative way the mappings involved in this work. (The detailed meaning of the symbols 
are given in the text with the specified equations.) 

(a) Continuum Hilbert space operator A onto phase space (Wigner function). Here a "line" underpins an Mutual 
Unbiased Basis (MUB) state projector \x' ,9){e.x'\. Eq.®. 

(b) Finite dimension, d=prime ^ 2 Hilbert space operator An onto points and lines of Dual Affine Plane Geometry 
(DAPG). Here a "point" a = m{b), underpins an MUB state projector \m,b){b,m\. 

(c) Finite dimension, d=primc 7^ 2 Hilbert space operator A^ onto points and lines of Affine Plane Geometry (APG). 
Here a "line", A = (m(— 1), to(0)), underpins A^IUB state projector \m,b){b,m\: 

il,)WAiq,p)^ Jdye^Py{q^y/2\A\q + y/2). Eq.® 

(b) y^„(j) = trANPf, j - a line. Eq.®. 

(c) V^^(e,77) = tr-(ijvQ5,„); (^,7?) - a point. Eq.^. 

Likewise we give now the Radon transform of the density operator mappings. 
ia.)tr{p\x',e){0,x'\) = {nWp){x',e) = 

=J ^S{x' - qcosO -psine)Wp{q,p). 6{x' ~ qcosB ~ psinO) = W\.^,e)e,x'\{q,v)- Eq.l^. 
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(b)trpA„ = {nuVp){a = m{h)) = \ Y.j {trpP^)K^^f, A„,, = VaAj)- Eq.m- 
{c)trpBx - (7^^Vp)(A = 77(0) = 1 (i'^pQe.r,) A(5,,,),a; A(^,„),;, = Vb, (e, 77)- i??-®. 

Appendix A: The lambda function, Af^^j) 

To prove the alternative meaning of the Lambda function, Eq. (j6ip . consider 

tri„P, = i^ir(P,P,v) = A,-„. 

Here we used the expression for in terms of Pj, Eq. ip^ . Evaluating this gives. 



Where we used Eas. (|571 



tr{AaPj} = Acj, A,,, = r'^. .QED. (61) 

0, J 3 a. 



Appendix B: Fluctuation Distillation Formula 



Given, Eq^, Pj = J^ce-, - I and, Eq(l37|, p/ = I, implies 



Thus, 



Recalling that, Eq([24l), = Aa allows 



QED 



Appendix C: The relation P? ^/ ^(0% = / 



- i=m(-l),m(0) 

Recalling Eqns.dSOl 1311 IMl) we have 

U otherwise 

Squaring this matrix we have trivially 1 at n=n', n=c/2. The only non nil element of the n-th row [n ^ c/2) of the 
matrix is along the column n'=c-n, where it is given by w^" The only column having non nil element at the 

row n=c-n' is the n-th column with the element cj^"-"" )'"(o). Thus {Pf^^/2 m(o))"."' ~ QED. 
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